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ABSTRACT 



A generalization of exterior calculus is considered by allowing the partial derivatives in the 
exterior derivative to assume fractional orders. That is, a fractional exterior derivative is 
defined. This is found to generate new vector spaces of finite and infinite dimension, 
fractional differential form spaces. The definitions of closed and exact forms are extended 
to the new fractional form spaces with closure and integrability conditions worked out for a 
special case. Coordinate transformation rules are also computed. The transformation rules 
are different from those of the standard exterior calculus due to the properties of the 
fractional derivative. The metric for the fractional form spaces is given, based on the 
coordinate transformation rules. All results are found to reduce to those of standard 
exterior calculus when the order of the coordinate differentials is set to one. 
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I. INTRODUCTION 



In recent years exterior calculus has been generalized by basing it on various graded 
algebras, see for example Refs. 1 and 2. Other attempts at generalization are based on non- 
associative geometries, see for example Refs. 3 and 4. In this paper another attempt at 
generalization is made using fractional derivatives in the definition of the exterior derivative. 
That is, a fractional exterior derivative is defined. Having a fractional exterior derivative 
gives rise to the notion of coordinate differentials of fractional order. These in turn can be 
used to define vector spaces of fractional differential forms. This formalism is found to 
produce an infinite number of finite and infinite dimensional vector spaces associated with 

each point P e E" (n dimensional Euclidean space). 

In sections 2 and 3 a brief review of exterior calculus and fractional calculus will be given to 
fix notation and provide a convenient reference. They are by no means complete, but are 
sufficient for the purposes of this paper. Section 4 defines fractional form spaces based on 
the fractional exterior derivative. Basis sets are given for the new vector spaces and notation 
is fixed. In section 5 the definitions of closed and exact are expanded to include the 
fractional form case. Once defined, the notions of closed and exact forms are examined for 
these new vector spaces. In both cases the results reduce to those found in standard exterior 
calculus when the order of the coordinate differentials is set equal to one. Coordinate 
transformations are worked out for the fractional form spaces in section 6. The 
transformation rules are somewhat more complicated than for standard exterior calculus. 
They do, however, reduce to the usual transformation rules when the order or the coordinate 
differentials is set to one. Having found the coordinate transformation rule, a metric for the 
fractional form spaces is constructed. Metrical properties of these new vector spaces will be 
investigated in a later paper. 
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The convention in the literature is that the coordinate index is a super-script. For the topics 
presented in this paper it is more convenient for the coordinate index to be a sub-script 
rather than the traditional super-script. To avoid confusion with this, the summation 
convention will not be used in this paper. 

II. BRIEF REVIEW OF DIFFERENTIAL FORMS 

The calculus of differential forms is an elegant branch of pure mathematics and a powerful 
tool in applied mathematics. A clear introduction to the field, with emphasis on applications, 

is given in Flanders 5) . Vector spaces at a point P e E" (n dimensional Euclidean space) 
can be constructed out of expressions of the following type 



c a = Y adx, , ... 
one forms, ~i (1) 



B = > b„dx: a dx : , 
two forms, H ifc " 1 (2) 



'n' forms, CO = w dx t a dx 2 a L Adx n . (3) 

Where the are the Cartesian coordinates of E". The above sums are taken over all 
possible values of the indices with the constraint that 



dx l a dxj = -dxj a dx t . (4) 
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The functions a t , b ip etc. depend only on P and may be real or complex depending on the 
application. If a k-form, y , is multiplying an m-form, fi , the following would be true, 



yA/i = (-l)" B /iA y. (5) 

The result would be zero if k + m > n. The exterior product, a , is distributive, associative, 
and anti- symmetric. The dimension of the vector space of k-forms over P e E" is 



ykj 



n 



which is zero if k > n. For the purposes of this paper let F{k,k,n) denote 



k\(n-k)\ 

the vector space of k-forms over P e E" . The apparently redundant 'k' in the above 
notation will be needed later for the fractional form case, as there is some additional 
freedom. 

The exterior derivative is defined as, 



d = t dx ^4r- ( 6 ) 



The exterior derivative maps k-forms into k+1 -forms and has the following algebraic 
properties. Let /and X be k-forms, and fi be an m-form, then 

d(y + A) = dy + dk, (7) 

d(y A(i) = (dy)A(i + {-\) k y a d/n, (8) 



d(dy) = 0. (9) 
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The last identity is called the Poincare lemma. A form, 7 , is called closed if dy = . A 
form, 7 , is called exact if there exists a form, fi , such that dfl = J . The order of fi is one 
less than the order of 7 . Exact forms are always closed. Closed forms are not always 

exact. The interested reader should consult Flanders 5) or Lovelock and Rund 6) for further 
details. 

III. BRIEF OVERVIEW OF FRACTIONAL CALCULUS 

There are many books that develop fractional calculus and the various definitions of 
fractional integration and differentiation. The reader should consult Refs. 7, 8, and 9 for 
further details and applications. For the purposes of this paper the Riemann-Liouville 
definition of fractional integration and differentiation will be used. T(q) is the gamma 
function (generalized factorial) of the parameter 'q' (i.e. T(n + 1) = n\ for all whole 
numbers, 'n'). 



d q f{x) = 1 r Am 

(d(x-a)Y r(-*)J-(x-Q 



W+l 



Re(q) < 0, 



(10) 



d q f{x) 



d n 



{d{x-a)) v dx n 



r 

J CI 



\q—n+l 



T{n-q) Sa (x-£) 



Re(<?) > 
n > q (n is whole) 



(ID 



The parameter 'q' is the order of the integral or derivative and is allowed to be complex. 
Positive real values of 'q' represent derivatives and negative real values represent integrals. 



Equation (10) is a fractional integral and equation (1 1) is a fractional derivative. In this 
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paper only real and positive values of 'q' will be considered. Notice that the derivative 
written in this form becomes a non-local object. Fractional derivatives have many 
interesting properties. For example the derivative of a constant need not be zero (the initial 
point 'a' in the above definition is set to zero in the following). 



d q l x~ q 



{dx) q Y{\-q) 



(12) 



The derivative of powers of x is 



dv _ r(p + i) P >-i 



(dxf T{p-q + \) q >0 



(13) 



Composition and product rales for fractional derivatives are given below. In the following, 
'n' is a whole number and 'q' is a complex number whose real part is greater than zero. 



d" d q d n+q ( . 

f\ x ) = 7^7 ^Jht/W 



dx n (d(x-a)) q {d(x-a))" 



— ^—^M = f(x) 



(14) 



(d(x-a)) q (d(x-a))~ q (15) 



d~ q d q 



(d(x-a))~ 9 (d(x-a)) 9 ^ 
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Composing derivatives where both have fractional order is given by the following formula. 



2q 2p+q k -\q-j 

' -/(*) = ,,, -E , 3 ,„ 



[d{x - a)) 1 ' [d{x - a)) q (d(x-a)) pci j=1 (d(x-a)) 



(x-a)- p - j 

r(i-P-j) 

(17) 



Where, 0<k-l<q<k , p > 0, and k is a whole number. The product rule is, 



W if8) = ^y. 



K dx j J 



(18) 



The above formula and definitions can be found in Refs. 7, 8, and 9. 
IV. FRACTIONAL FORM SPACES 

If the partial derivatives in the definition of the exterior derivative are allowed to assume 
fractional orders, a fractional exterior derivative can be defined. 



d v = J^dxJ 



!=1 



(d(x t -a)) V 



(19) 



Note that the subscript 'i' denotes the coordinate number, the superscript ' v' denotes the 
order of the fractional coordinate differential, and a i is the initial point of the derivative. 



Sometimes the notation d v , will be used to denote 



(d(x i -a,)) v 



7 



In two dimensions (x,y), the fractional exterior derivative of order v of x p , with the initial 
point taken to be the origin, is given by 



d v x p = dx v r ( p + 1 ) r p ~ v + dv v — no\ 

T(p-V + l) y yT(l-v)' (20) 



For specific values of the derivative parameter the following results are obtained. 

v = d°x p = 2x p (21) 

v = l d l x p = dx'px 1 " 1 ( 22) 

V = 2 d 2 x p = dx 2 p(p - l)x p ~ 2 (23) 

By analogy with standard exterior calculus, vector spaces can be constructed using the dxj . 

Let F( v ,m,n) be a vector space at P e E" • 'V denotes the sum of the fractional 
differential orders of the basis elements, 'm' denotes the number of coordinate differentials 
appearing in the basis elements, 'n' the number of coordinates, and {x t } are the Cartesian 

coordinates for E". For example, a basis set for F(v,l,n) would be \dx\, dx v 2 , . . . , dx v n } 
and an arbitrary element of F(v,l,n) would be expressed as 



a = Y,a,dxr. (24) 
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For a fixed v this is an 'n' dimensional vector space. Also note that there is a different 
vector space for each value of v. For V = 1 the one forms from exterior calculus are 
recovered. Now suppose that the basis elements are made up of two coordinate 
differentials, F(v,2,n) . In this case the basis set is more complicated. 



{dx? n Adx? 21 , dxf" Adx? 1 , . . . , dx^~ lm Adx^ \fi y + /l kj = v} 



(25) 



Note that dxf" a dx^ 2i would be zero if and only if fi n = fi 2l ,etc. An arbitrary element of 
F(v,2,n) would be expressed as a sum of the form 



p=£ I I v (a 7 (v„v-v ( X' a^p) dv, 



(26) 



where /3 (! (jU,ju) = 0. Unlike the previous vector space, F(v,l,n), this is clearly infinite 
dimensional for any value of v. Not only is it infinite but it is uncountably infinite. An 
arbitrary element of F(v,3,n) would be expressed as an integral of the form 



n n n 



P = IIL\ V l V " V '(Av,(v,v-v.,v-v ( -v.X' Adx]> Adxr'-^dvjdv,. 



(27) 



With each step up on the middle index of F(v,m,n) another integral and summation is 
included. A basis set for this vector space would be 



F(v,m,n) = \dxp' a dx- i2> a ... a dx-"" 1 , . . . 



XX- 



(28) 
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The basis elements range over all possible combinations of the fractional coordinate 
differentials and all possible choices for the fjC s. Note that 'm' need not be less than or 
equal to 'n' . 

Let P e E" , and let A e F(v,m,n) and B e F{ji,k,n) at the point P. Then the exterior 
product of A and B maintains the antisymmetry property of equation (5), 

AaB = {-l) km BAA e F(fi + v,k + m,n). (29) 

If k + m> n, AaB need not be zero. Equation (7) is also maintained due to the linearity 
of the fractional derivative. Equation (8) is not maintained due to the product rule for the 

fractional derivative (see equation (18)). Note also that d v maps F{fi,k,n) into 

F{fi + v,k + l,n) . 

V. CLOSED AND EXACT FRACTIONAL FORMS 

By analogy with exterior calculus the notions of closed and exact can be extended to 
fractional forms. 

Let g g F(n,k,n) then g is v -exact if 3 an / e F{ji-v,k-\,n) such that d v f = g. 
Let g e F(v,k,n) then g is H -closed if d M g = 0. 

To examine the notion of exactness (integrability conditions) in the fractional form case the 
kernel is needed for the fractional derivative operator. This will be denoted by, Ker[d^ . In 
the following the initial point of the derivative will be taken to be the origin. Solve 
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dj(h) = 0. 



(30) 



Equation (30) is solved using equation (1 1). Let 'm' be the first whole number greater than 
or equal to ' v ' then 



m— 1 

\V-m / \k 



h = (x) v ~ m £q (*,■) 



(31) 



k=0 



This is basically the result from Oldham 7) page 155. The c k ' s can be functions of the 
other coordinates. The kernel for the operator, d v (when restricted to act only on scalar 
functions) is similarly constructed, Ker(d v } . 



m ~)V r 

d v f= Ydxj^- = o 



(32) 



( n \ 



f= n^, 

V 1=1 J 



m—l m — 1 



(33) 



The C k L k are now constants and 'm' is once again the first whole number greater than or 
equal to ' v ' . 

The fractional integrability conditions can now be constructed for the following restricted 
case. Let g be a v - form in F(v,l,n) 



i=i 



(34) 
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When can a - form, f, be found such that d v f = g ? If such an f exists it will be contained 
in the family of functions given by 



m— 1 



(35) 



k=0 



Recall that there is no sum over the repeated indices, except where a summation symbol is 
encountered. This solution must satisfy, djf = a 7 , 



m-l \ 



*=0 / 



(36) 



The above equation is to be solved to determine the unknown functions c k , and must be true 
for all values of 'i'. Equation (36) can be rearranged to give 



m-i 



k=0 



(37) 



The left hand side of (37) is a polynomial of order 'm-l' in the variable x n hence it can 
only be solved for the c k if the following is true 



a, 



dx" 



= 



(38) 
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The above equation is the integrability condition for fractional forms of type F(v,l,n) • 

Note that if v = m = 1 the usual integrability conditions from exterior calculus are 
recovered. 



Consider the second definition and examine what it takes to be closed for fractional forms. 

n 

Let a ~ L* a i^ x i G ^^ V '^' n ^ and consider its fractional exterior derivative. 



d»a = £ d"(a,. dxj) 



i=\ 



(39) 



i=\ 7=1 k=0 



k 



V dx J 



(40) 



In the last sum of equation (40) 'k' takes on only whole number values hence 



dx) 



(dx]) = V k > 1. 



(41) 



This reduces equation (40) to the following 



i=l 7=1 



a. 



(42) 



Since and dx] are linearly independent, provided fi ^ v or z ^ j, J M a = if and only 



if 
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j-YCC, = 0. (43) 



In other words a i e Kerid 1 } ). For the special case of ji = v the symmetry from equation 
(29) can be used to obtain 



For v = 1 the usual result from exterior calculus is recovered. 



^-aj - l-a,. = (45) 
OX; ox , 



VI. TRANSITION TO CURVILINEAR COORDINATES 

When coordinate transformation rules are worked out for exterior or tensor calculus (see 
for example Refs. 5 and 6) the following construction can be used. Let {x { } and be 
two coordinate systems with a one to one mapping between them in some neighborhood of 
P e E" . Take to again be Cartesian coordinates and {y t } to be curvilinear 

coordinates. Assume the {x { } can be written smoothly in terms of the {y t } . 



x 



i =x i {y) (46) 



The exterior derivative is then applied to the above expression giving the following. 
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dx i = dx,{y) (47) 



X>*|r = !><!r (48) 

7~1 dx k dy, 



£>,<5f = £dy,fL (49) 
ti Tt dyi 



dx t = Vdy^ (50) 



This is very straightforward and can be adapted to the fractional form case. In the two 
coordinate systems the fractional exterior derivative d v takes the following forms 



d v = idx;—^—f, (5D 

,-=i \d(x i -a,)) 



and 



,-=i (o>(v, -a,)) 



Where a, is the initial point of the derivative in the Cartesian system and a i is the same 
point but in the curvilinear coordinates. Recall that x i is in the kernel for the operator 
d I dx k for k # i, and when k = i the result is one. For the fractional case the same type of 
object is needed. Consider a function a k that maps points in E" into the complex numbers 

15 



a, = 



r(D 
r(v+i) 



w* J 



(53) 



a i is the initial point for the fractional derivative. The function a k was chosen so that it 



would be in the kernel of 



(<?(*,. - a,)) 1 



for i ^ k, and for i = k 



d v a„ 



(d(x k -a k )) V 



1 . 



(54) 



If the fractional exterior derivative is applied to a k in the two different coordinate systems 
the following coordinate transformation rule can be obtained (where all quantities on the 
right hand side of (55) must be expressed in terms of the {y t } coordinates). 



dx v k 



dy? d v 



trr(v + i)(^-5,.)) 1 



!=1 



(55) 



Note that for v = m = 1 the usual coordinate transformation rule is recovered. The 
coordinate transformation matrix for the fractional forms will be denoted by 



J-(x,y,v) = 



nv + D(d{y-~a)) V 



f / \V-m 

U( x j~ a j) 



(56) 



dxl = £dyJJ*(x,y,v) 

7 = 1 



(57) 
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A fractional form, 



A = i,A k (x)dx v k (58) 

k=\ 



written in the {*, } coordinates would be transformed into the {y t } coordinates according to 



A = LL A k {x(y))j*(x,y,v)dy! . (59) 

k=l i=l 



Reversing the coordinate transformation yields 



A = LLL Mx)JXx>y>v)Jj(y,x,v)dx] , (60) 

k=\ i=l j=l 



^>5*=£/;(y,^v)/*U,y,v). (61) 

As an example consider the coordinate transformation for two dimensional Cartesian 
coordinates to polar coordinates. The initial point for the fractional derivatives is taken to be 
the origin. 



jCj = rcos(6>) and x 2 = rsin(6>) (62) 



The coordinate transformations for the fractional differentials are then 
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r(2v-m + l) cos v (0) v 



, 2v-m ~\v 



v—m j v , 

r dr + 



' cos v (0) A 



r(v + i)r(v-m + i) sin m - v (0) rxv+i^y^skr-^J 



d0 v , (63) 



dxl = 



r(2v-m + l) sin v (0) v 



2v-m ")V 



r dr + 



r(v + l)r(v-m + l)cos m " v (0) T(v + l)(d0) 



yCOS m - V (d)j 



dG v . (64) 



For v = m = 1 the transformation equations from exterior calculus are recovered. 



dx l = cos(0)dr - rsm(9)d0 



(65) 



dx 2 = s'm(9)dr + rcos(0)dO 



(66) 



Having found the coordinate transformation rule a metric for F{y,l,n) can be constructed 
just as is done in exterior calculus (see page 46 of Lovelock and Rund 6) ) 



8ij(y> v ) = £/*U,y,v)/ 7 *U,y,v) 



k=l 



(67) 



Which can be used to give a fractional line element 



ds2v = t.Sijiy^dy] ® dy] 



(68) 



or, 



dsV = \t,s tJ (y,y)dy] ® dy] 



(69) 
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Where ® is the symmetric product for coordinate differentials. 



VII. CONCLUSION 



In this paper a natural extension of the exterior derivative to a fractional exterior derivative 
was considered. It was found to generate some new vector spaces, both finite and infinite in 

dimension. Of particular interest is the observation that at each point P e E" there are an 
infinite number of n dimensional vector spaces. Due to their similarity with tangent spaces, 
perhaps an appropriate name for these new vector spaces should be fractional tangent 
spaces. Notions of closed and exact were also defined for these fractional form spaces. 
Integrability and closure conditions were investigated for the special case of F{y,l,n). The 
results produced were found to reduce to the standard results from exterior calculus when 
the order of the fractional exterior derivative was set equal to one. Coordinate 
transformation rules were also worked out for F(v,l,n). The transformation rules are 
somewhat more complicated than for traditional exterior calculus. They do however reduce 
to the usual formula when the order of the coordinate differentials are set equal to one. 
Coordinate transformation rules give rise to a metric for F(v,l,n). Properties of the metric, 
such as its associated covariant derivative, will be investigated in a later paper. 
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